Fast entanglement detection for unknown states of two spatial qutrits 
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We investigate tiie practicality of the method proposed by Maciel et al. [Phys. Rev. A. 80, 032325 
(2009)] for detecting the entanglement of two spatial qutrits (3-dimensional quantum systems), which 
are encoded in the discrete transverse momentum of single photons transmitted through a multi-slit 
aperture. The method is based on the acquisition of partial information of the quantum state through 
projective measurements, and a data processing analysis done with semi-definite programs. This 
analysis relies on generating gradually an optimal entanglement witness operator, and numerical 
investigations have shown that it allows for the entanglement detection of unknown states with a 
cost much lower than full state tomography. 

PACS numbers: 03.67.Mn, 42.50.Dv 



I. INTRODUCTION 

One of the main goals of quantum information the- 
ory is to find general and simple procedures that deter- 
mine whether a composite system is entangled or not [1]. 
The simplest example is the entanglement detection of 
pure bipartite states of qubits (2-dimcnsional quantum 
states). In this case, the entanglement of the compound 
state can be determined directly from the reduced den- 
sity matrixes [2, 3]. This approach can be generalized for 
distinct types of quantum states [4-6], and they all rely 
on state dependent nonlinear functions that character- 
ize entanglement quantitatively. These functions can be 
evaluated once one has the knowledge of the compound 
density operator. 

Unfortunately, the experimental determination of com- 
posite states can require an extremely fastidious work 
[7, 8]. An unknown D-dimensional quantum state (qu- 
dit) is represented by a D x density operator, which 
requires — 1 independent real parameters for its spec- 
ification. The state reconstruction is done with the tech- 
nique of quantum tomography (QT), where these param- 
eters are associated, usually, with an overcomplete set 
of observables that are evaluated on single copies of the 
quantum system in consideration [9, 10]. Besides of the 
large amount of measurements necessary in the QT, it 
may also require sophisticated programs for doing the 
unavoidable numerical optimization of the acquired data 
[11]. Depending on the dimension of the composite sys- 
tem, this optimization may be even impossible with to- 
day's technology [11]. 

Therefore, it is legitimate to ask if there are more 
practical schemes that allow for the entanglement de- 
tection of composite states, which do not require the 
full state reconstruction. One possibility is to consider 



'Electronic address: KlimaCSudec.cl 



the so-called collective measurements [12-18], where the 
idea is to detect the entanglement by measuring observ- 
ables that act on simultaneously available copies of the 
composite state. However, the simultaneous generation 
of copies of a quantum state may not be possible in 
some circumstances, and another solution for the entan- 
glement detection is the use of entanglement witnesses 
operators [19, 20]. These operators are observables {W) 
with non- negative mean values for all separable states, 
such that the measurement of a negative mean value 
[rr[W pab) < 0] witness the entanglement of the com- 
pound state (pab)- It is also worth to note that, recently, 
a link between collective measurements and single-copy- 
entanglement-witnesses has been presented [21]. 

In general, there is no universal witness operator (see, 
however, Ref [22] for a particular case) and one must 
assume a prior knowledge of the composite state for con- 
structing the optimal EW operator [23, 24]. Nevertheless 
this approach is of limited application, since this prior 
knowledge of the state may not be available. Moreover, 
the witness operator will give reliable information only if 
the considered state is detected. Otherwise one can not 
distinguish whether the state was separable or contained 
a distinct type of entanglement. 

The efficient construction of the entanglement wit- 
nesses is also another important drawback for this tech- 
nique of entanglement detection. There are many dis- 
tinct techniques presented [25-29], but the optimal way 
for witnessing the entanglement of a given state is still 
an unsolved problem. 

For circumventing the practical problem of witness- 
ing the entanglement of unknown states, there has been, 
recently, the development of a numerical and interac- 
tive method, where the witness operator is constructed 
and evaluated while the experimental data is collected 
[30, 31]. For detecting the entanglement of a compound 
state from an incomplete set of measurements, one can 
also use the techniques of [32, 33]. The need for such 
procedures arises often in real applications, and in this 
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work we present the first experimental demonstration of 
the method described in [30, 31]. It is interesting to note 
that numerical investigations have shown that it allows 
for the entanglement detection with a cost much lower 
than QT. 

For testing the method of Maciel et al. [30, 31], we gen- 
erated bipartite entangled states of 3-dimcnsional quan- 
tum systems (qutrits). The qutrits are encoded in the 
discrete transverse momentum of single photons trans- 
mitted through a multi-slit aperture, and we refer to 
them as spatial qutrits [34]. Pure entangled states with 
distinct degrees of entanglement [2, 3] were used to study 
the method when it is applied to a system with a variable 
resource, namely entanglement. We also considered dif- 
ferent measurement settings for each state generated, and 
this allows one to verify the experimental convergence of 
the method under distinct scenarios. 

The experimental test of a theoretical method is im- 
portant to determine the practicality and limitations of 
the method in actual implementations. Besides, it is im- 
portant to note that the work presented here can also be 
interpreted as a new technique for estimating the entan- 
glement of transverse spatially correlated quantum sys- 
tems, whose characterization has been the theme of re- 
cent investigations [35-39]. As it was shown in [33, 40], it 
is possible to associate lower bounds for some entangle- 
ment measures, given that the mean value of a witness 
operator has been determined. 



II. EXPERIMENT 

A. Generating the spatial qutrits 

In Fig. 1, it is represented the setup used for the gen- 
eration and entanglement detection of the spatial qutrits 
states. A pigtailed single mode diode laser, operating 
at 670 nm, is attenuated to the single photon level with 
absorptive neutral density filters (NDF). The particle be- 
havior of the light was confirmed by using an avalanche 
photodetector (APD) to register the transmitted pho- 
tons. After 2 m of free propagation (which is the over- 
all extension of our experimental setup), the maximal 
single count rate recorded was 50000 s~^. Taking into 
account the coherence length of the diode laser and the 
APD properties, this shows that, on average, fewer than 
10"'^ photons is presented in the experiment at any given 
moment. 

After this initial test, the transmitted photons are sent 
to a transmissive spatial light modulator working for 
amplitude-only modulation. It is composed of two po- 
larizers (P) and a twisted nematic liquid crystal display 
(LCD). The polarizers are placed close to the LCD panel 
(0.5 cm away from each side). The LCD is connected 
to a PCB interface which allows its control from a com- 
puter. This SLM is addressed with a 9-multi slit. The 
slit's width is 2a = 104 /xm and the distance between 
two consecutive slits is d = 208 /im. In this case, the 
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FIG. 1: Experimental Setup. See the main text for details. 



overall SLM transmission is 20%. The state of the single 
photons transmitted through this SLM can be written 
as [34, 41] 



(1) 



where the state \l) is defined by \l) = 
Jdqe~'^'''-'^smc{qa)\lq). It represents the state 
of the photon transmitted by the /th-slit of the SLM. 
\lq) is the Fock state of a photon with the transverse 
wave vector q. These states 1 1 ) form an orthonormal 
basis in the 9-dimensional Hilbert space of the trans- 
mitted photons and, therefore, they are used to define 
the logical spatial states. The state of Eq. (1) represents 
a 9-dimensional quantum state that is encoded in the 
discrete transverse propagation modes of the transmit- 
ted photons. Alternatively, it can also be thought as a 
composite system of two qutrits, which can be entangled 
or not. For seeing this, one just need to consider a new 
label for the |/) states, like for example |/) — )• 
such that: |-4) -> jOO); |-3) ^ jOl); |-2) |02); 
1-1) ^ jlO); |0) ^ 111); |1) ^ |12); |2) ^ |20); 
I 3) ^ I 21) and I 4) ^ I 22). In this case, the state of 
the transmitted photons is given by 



2 

i,j=0 



(2) 



These spatial qutrits are encoded in the transverse spatial 
modes of a single particle, and therefore they can not be 
used to test quantum non-locality [42, 43]. Nevertheless, 
they can be used to study other fundamental aspects of 
quantum mechanics [44]. 

To generate spatial qutrits states with different degrees 
of entanglement, it is necessary to modify the coefficients 
aij of Eq. (2). These coefficients are dependent on the 
spatial profile and the wavefront curvature of the laser 
beam at the SLM plane. They are also dependent on the 
slit's transmissions addressed in the SLM, and as it was 
shown in [41], the spatial light modulator can be used to 
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modify, independently, each slit transmission. This mod- 
ulation is done by changing the grey level of the LCD 
at the points where the slits are addressed, and thus, 
it allows for a controlled generation of distinct types of 
spatial qutrits entangled states. This manipulation has 
been demonstrated to be coherent, such that the modi- 
fied states can still be described by a coherent superpo- 
sition [41]. 

In our work, the qutrits states were generated by send- 
ing a well coUimated beam to the first spatial light mod- 
ulator. In this case, the coefficients aij are expected to 
be real, and the initial state of Eq. (2) can be modified 
by the SLM to 



odl 



ij=0 



(3) 



where = y/^, with tij representing the transmission 
of the ij-slit. is a normalization constant. In Fig. 2(a) 
it is shown a comparison between the laser spatial profile 
at the plane of the SLM and the 9-multi slit being ad- 
dressed on it. This graph shows how each slit of the array 
was initially illuminated. In Fig. 2(b) one can see the de- 
pendence of the SLM transmission and the grey level of 
its LCD. In Fig. 2(c), Fig. 2(e) and Fig. 2(g) it is shown 
the modulations used for generating three distinct types 
of spatial qutrits states. The corresponding expected am- 
plitudes of these states are shown in Fig. 2(d), Fig. 2(f) 
and Fig. 2(h), respectively. These amplitudes are calcu- 
lated by taking into account the spatial laser profile and 
the transmissions of the slits being addressed in the SLM. 
As it will be discussed on next section, these states have 
different degrees of entanglement. The first state gen- 
erated is a highly entangled state of two spatial qutrits, 
while the second and the third states are generated by 
reducing gradually the degree of entanglement between 
the qutrits. 
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FIG. 2: The expected spatial qutrits states. In (a) it is shown 
a comparison between the laser beam spatial profile and the 
9-multi slit aperture addressed in the first SLM. In (b) it 
is shown the dependence of the SLM transmission and the 
grey level of its LCD. In Fig. 2(c), Fig. 2(e) and Fig. 2(g) it 
is shown the modulations used for generating three distinct 
types of spatial qutrits states. The corresponding expected 
amplitudes of these states are shown in Fig. 2(d), Fig. 2(f) 
and Fig. 2(h), respectively. 



B. Characterizing the states generated 

Even though we can estimate which are the spatial 
qutrits entangled states generated in our experiment (See 
Fig. 2), this is not enough for demonstrating that the 
method of [30, 31] allows for a fast detection of their 
entanglement. For doing such study it is first necessary 
to completely determine these states, so that we know 
precisely what we should obtain in the test. 

In this work we use the recently proposed technique 
of Ref [45] for performing the quantum tomography of 
the spatial qutrits states. It is based on projecting these 
states onto the vectors of the mutually unbiased bases 
(MUBs), and it has the advantage that the overcomplete 
set of vectors used for the reconstruction has a number 
of projections that is minimal [46, 47]. In the MUBs, any 
two vectors of different bases have the same overlap's ab- 
solute value. There exist 10 MUBs for a 9-dimcnsional 



Hilbert space. The existence of these MUBs and con- 
structive routines for obtaining them, has been given by 
Wootters in [47]. The density operators of the spatial 
qutrits states can be represented in terms of the MUBs 
as [46] 

/^^EEp^^'nL")-/, (4) 

Q— 1 m— 1 

where vt^ = \it^)i(il^t\. and vt^ = rr(pn,^;?^) is the 
probability for projecting p onto the m-th state jV'm^) of 
the a MUB. We note that the states jV"™ are given ex- 
plicitly in Ref [47]. / represents the 9-dimensional unity 
operator. Thus, the QT based on the MUBs of the spa- 
tial qutrits is implemented by measuring the probabilities 
Pm' of the 90 corresponding projectors (Ilm^)- 

For projecting the generated spatial qutrits states onto 
the vectors of the MUBs, it is necessary to consider a sec- 
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FIG. 3: Reconstructed states. In (a), (d) and (g) it is shown the recorded probabilities and the expected ones (insets) that are 
calculated from the states given in Fig. 2(d), Fig. 2(f) and Fig. 2(h), respectively. In (b) [(c)], (e) [(f)] and (h) [(i)] it is shown 
the real (imaginary) parts of the corresponding reconstructed density operators in the logical base. On the insets it is shown 
the expected parts of these density operators. 



ond SLM and a spatial filtering process as it is described 
in details in [45] . Briefly, a second SLM is set with its liq- 
uid crystal display (LCD2) placed at the plane of image 
formation of the first LCD. In our case we used a 7.5 cm 
focal length lens (LI) (See Fig. 1). This SLM must be 
configured for doing phase-only modulation and then it is 
used to modulate the phase of the transmitted light, in- 
dependently, at each point of slit image formation. This 
is done by changing the gray level of the LCD2, indepen- 
dently, at these points. At the plane of image formation, 
the diffracted photon is again described by Eq (3) [48] 
and, therefore, the modulation being done by the sec- 
ond SLM allows for the modification of the imaginary 
parts of the spatial qutrits states. After the reflection on 
this SLM, the light is coUimated and focused by a set of 
lenses (L2 and L3) at the detection plane, where a point- 
like detector (APD) is used to record the single photons 
at the center of the interference patterns formed. This 
detector is composed of a conventional avalanche photo- 
counting module with a slit of 20 fj,m in front of it. The 
spatial filtering being introduced by the set of lenses and 
the point-like detector, together with the phase modula- 



tions of the second SLM, allows for projecting the spa- 
tial qutrits states onto the |'0m ■*) MUBs vectors of their 
9-dimensional Hilbert space [45]. 

The recorded probabilities Pm\ for the first, second 
and third state generated in our experiment are shown in 
Fig. 3(a), Fig. 3(d) and Fig. 3(g), respectively. On the in- 
sets of these figures it is shown the expected probabilities, 
which are calculated from the states given in Fig. 2(d), 
Fig. 2(f) and Fig. 2(h), respectively. In Fig. 3(b) and 
Fig. 3(c) it is shown the real and the imaginary parts of 
the reconstructed density operator of the first state. On 
the insets it is shown the expected real and imaginary 
parts. For obtaining this density operator we replace 
the values of the probabilities p™-* shown in Fig. 3(a) 
into Eq. 4, and then optimize the resulting matrix to 
guarantee that it is positive semi-definite [10]. For do- 
ing this optimization we considered the "forced purity" 
approach, which has been shown to be sufficient when 
one is dealing with high-dimensional quantum systems 
of low-entropy [11]. As it was shown in [45], this is ex- 
actly our case. The "forced purity" numerical technique 
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generates a pure quantum state for which the mean val- 
ues of the measurements considered are very close to 
the ones obtained in the experiment. The fidelity [49] 
of the first reconstructed state with the expected one 
[Fig. 2(d)] was F,tatei = 0.84 ± 0.04. The real and the 
imaginary parts of the density operator obtained for the 
second (third) state arc shown in Fig. 3(e) [Fig. 3(h)] and 
Fig. 3(f) [Fig. 3(i)], respectively. The expected parts of 
these density operators are also shown on the insets of 
these figures. For the second state, the fidelity between 
the obtained density operator and the expected one was 
Fstate2 = 0.87 ±0.04. For the third state we got a fidelity 
of Fstates = 0.92 ± 0.04. The high value of the fidelities 
between the reconstructed states and the expected ones 
is usually accepted as a demonstration that the experi- 
mental setup is working as it is intended to, so that we 
may use it for further studies. 

Next we investigate the numerical technique of [30, 31]. 
This technique is dependent of the degree of entangle- 
ment of the composite system and so, it is important 
that we determine also the degree of entanglement of the 
reconstructed states of our experiment. For calculating 
the degree of entanglement of the spatial qutrits states 
we use the / concurrence introduced by Rungta et al. [6] , 
which is a generalization of the concurrence [3] to arbi- 
trary dimensional bipartite composite systems. For pure 
bipartite systems, it is dependent only on the purity of 
the marginal density operators 

(Ci"i)2 = 2[l-MPD], (5) 

where k — A^B denotes the subsystems. In this case 
of pure bipartite systems, the / concurrence covers the 
range of [50] 

< (Ci"i)^ < ^i^, (6) 

where n = min[nA,nB], with ua and ub defining the 
subsystems dimensions. In the case of pure spatial qutrits 

states we have that < (C^^)^ ^ |- ^'^^^ Eq. (5) we can 
calculate the / concurrence of the reconstructed spatial 
qutrits states given in Fig. 3. The calculated values are 
shown in Tab. I. 



/ concurrence 



State 




1 


1.20 


2 


0.45 


3 


0.27 



TABLE I: / concurrence of the reconstructed spatial qutrits 
states. 



C. Fast entanglement detection for unknown 
spatial qutrits states 

Overview of the method — The method of Macicl 
et al. [30, 31], can be seen as a variational approach that 
allows one to estimate some physical properties of an 
unknown quantum system. It requires the acquisition of 
partial information of the state of the system, obtained 
through projective measurements, and it consists of a 
interactive routine, where semi-definite programs [51] are 
used to gradually estimate the expectation value of an 
observable while the data is collected. 

This estimation has been showed to be optimal and 
reliable [31]. Optimal in the sense that it converges to 
the correct mean value of the observable in considera- 
tion. Reliable, because the program will not attribute to 
a system an upper value for its physical property. The 
obtained expectation value will always be smaller in ab- 
solute value than the true one. This is an important 
property, specially, when one is detecting the entangle- 
ment of the system with entanglement witnesses. In this 
case, the entanglement of the system will never be over- 
estimated. 

The basic idea of the method is the following. Suppose 
some projectors have been measured on the unknown 
quantum state. These projectors can be, for example, 
some of the projectors we used above (Ilm ^ ) , in the case 
one is dealing with spatial qutrits states. Once these pro- 
jectors have been measured, one will have determined the 
corresponding projection probabilities (some of the prob- 
abilities pm'' in our case), and it is then possible to write 
the expectation value of an observable as a function of 
two terms: one which is a function of the known prob- 
abilities, and one which depends on the unknown ones 
[31]. Thus, the idea is to perform a variational estima- 
tion of the expectation value, by doing an appropriate 
minimization of the term which depends on the unknown 
probabilities. This problem consist of a linear optimiza- 
tion problem, which can be solved efficiently with semi- 
definite programs. 

Since the method allows for the estimation of observ- 
ables, it can be used for doing the entanglement detection 
of composite states with the technique of entanglement 
witnesses operators [30]. To build the witness operator, 
the strategy used is the one of Ref [26] , which allows for 
the construction of an optimal trace one witness. Be- 
sides, the method can be used for the entanglement de- 
tection of unknown states. This is possible because the 
witness operator is built according to a "guess" state, p, 
which is efficiently inferred from the projective measure- 
ments already performed in the composite system [31]. 
When the number of projective measurements increases, 
we have that p — >■ p (where p represents the real/actual 
state of the system). When the measurements are equiv- 
alent to the full QT, the witness operator corresponds to 
the optimal one for the compound system. 



6 



(a) 



(b) 



MUB Order: 

tt=[10, 1, 2, 3, 4, 5, 6, 7. 8, 9] 



20 30 40 50 60 

Measmed MUB ProjectoK 



'Ci 0- 



MUB Order: 

a=[10, 1, 2, 3, 4, 5, 6. 7, 8, 9] 



10 20 30 40 50 60 70 80 90 

Measured MUB Projectors 



0.1 










MUB Order: 






1-0.1 




a=[10, 1, 2, 3, 4, 5, 6, 7, 8, 9] 










-0.2 





10 20 30 40 50 60 70 80 90 

Measured MUB Projectors 

FIG. 4: Entanglement detection for the generated spatial 
qutrits states, without assuming any knowledge of them. In 
(a) it is shown the results obtained for the first state gener- 
ated in our experiment. In (b) [(c)] it is shown the results 
obtained for the second (third) state generated. See the main 
text for details. 



Fast entanglement detection for unknown spa- 
tial qutrits states — Now we show how the program 
was used for doing the fast entanglement detection of 
the states generated in our experiment. No prior knowl- 
edge of these states is assumed. The data used in the 
analysis is the same acquired during the state recon- 
struction. The recorded probabilities for these states are 
shown in Fig. 3(a), Fig. 3(d) and Fig. 3(g). While the 
data was being collected, the program was used to detect 
the entanglement of the generated states. In Fig. 4(a), 
Fig. 4(b) and Fig. 4(c) it is shown the obtained results 



for the first [ C 



,(3) 



= i.2o; 

2 



second [(C^AB^ = 0.45] 



ing the convergence observed for each state. For the 
first state (Fig. 4(a)), the convergence was faster. 23 
measurements were necessary for a reasonable estimation 
of the witness expectation value. For the second state 
(Fig. 4(b)), which is less entangled, the system required 
around 30 measurements. For the last state (Fig. 4(b)), 
which has the lowest / concurrence, more than 55 mea- 
surements were necessary. This is important because as 
it was shown in [33, 40], it is possible to associate lower 
bounds for some entanglement measures, given that the 
mean value of a witness operator has been determined. 
Thus, we may say that the method of [30, 31] also al- 
lowed for the fast estimation of the entanglement of the 
spatial qutrits states. Moreover, the behavior observed 
is exactly what one should expected, since it should be, 
in principle, easier to estimate the entanglement of states 
that are more entangled. 

Further analysis — As it is discussed in details in 
[31], the numerical technique we are investigating turns 
out to be also an algorithm that allows for QT with an 
incomplete set of observables. As we mentioned above, 
this is done with a numerical estimation of the state of 
the system, that takes place after some projective mea- 
surements have been performed. In Fig. 5, it is shown 
for the first state of our experiment, a comparison be- 
tween the state generated by the method of Maciel et al. , 





22 00 



and the third state [{c^ab) — 0-27], respectively. The (c) 
order of the measurements performed is indexed with 
a = [10, 1, 2, 3, 4, 5, 6, 7, 8, 9]. This means that the first 9 i 
projective measurements implemented correspond to the o.b 
projectors of the lOth-MUB, the 10th to 18th measure- .^o.e 
ments correspond to the projections onto the Ist-MUB 1 0.4 
and etc... The lOtli-MUB, as one can see in Fig. 3, cor- 0.2 
responds to the logical base of the qutrits Hilbert space, 
and this is why it was measured first. This was done 
according to the procedure described in [45]. 

From Fig. 4 one can observe that the entanglement 
was always detected much faster than QT. While QT 
requires a total of 90 measurements to determine the en- 
tanglement of the generated states, the method of [30, 31] 
allowed the entanglement detection already on the sec- 
ond measured MUB. This happened independently of the 
degree of entanglement between the spatial qutrits. 
Another important observation to be done is regard- 
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FIG. 5: State reconstruction for the first entangled state of 
spatial qutrits. In (a) and (b) it is shown a comparison be- 
tween the final obtained state pi, and the expected one (in- 
sets) calculated from the amplitudes shown in Fig. 2(d). In (c) 
it is shown the fidelities between these states while the mea- 
surements were being performed. In (d) it is shown the fideli- 
ties between these states when the measurements are taken 
in a distinct order. 
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and its expected form calculated from the amplitudes 
given in Fig. 2(d). One can see that the reconstruction 
presented a high value of fidelity, and this can be seen 
as an evidence of the good quality of the reconstruction 
done. The same behavior has been observed for the other 
states. The fidelities between the final density operators 
obtained with this method, for the second and third state 
generated in the experiment, while compared with their 
expected form, were 78% and 82%, respectively. 

It is interesting to note in Fig. 5(c) and Fig. 5(d) that 
the convergence of the method depends on the measure- 
ment settings. In Fig. 5(c) it is shown the fidelities ob- 
tained while the measurements were being performed. 
Figure 5(d) was generated by applying the method to 
the data after it has been organized in a different order. 
One can see that the convergence was faster in Fig. 5(c). 

Since we have seen that the convergence of p changes 
with the order of the measurements performed, it is in- 
teresting to have a look on how the expectation values 
of the witnesses operators are affected. This is shown in 
Fig. 6. In Fig. 6(a) and Fig. 6(b) we have the entan- 
glement detection for the first state, when the measure- 
ments are considered in a different order at the program. 
In Fig. 6(c) and Fig. 6(d) we have the entanglement de- 
tection for the second state. In Fig. 6(e) and Fig. 6(f) we 
have the same type of analysis for the third state. One 
can see that even though the convergence of the method 
is lower than before (Fig. 4) , it can still detect the entan- 
glement of the generated states much faster than QT. 



III. CONCLUSIONS 

In this work we have presented a detailed investiga- 
tion of the technique recently introduced by Maciel et al. 
[30, 31], which allows for the fast entanglement detection 
of unknown states of composite systems. For testing the 
method we considered entangled states of two qutrits, 
which were generated with distinct degrees of entangle- 
ment. This allowed the experimental test of the method 
under distinct scenarios. The qutrits states were encoded 
on the discrete transverse modes of single photons trans- 
mitted through a diffractive aperture. Our results show 
that the method of [30, 31] can indeed provide fast entan- 
glement detection even when experimental errors, such 
as misalignment and imperfections, are present. Even 
though we have studied the behavior of the method us- 
ing only pure entangled states, it has been shown to be 



also valid for mixed states. It is also important to note 
that the method can be applied to any set of projective 
measurements used for doing the quantum tomography in 
previous experiments. The work presented here can also 
be seeing as a new technique for doing the entanglement 
estimation of spatially correlated quantum systems. 
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FIG. 6: Entanglement detection for the generated spatial 
qutrits states, without assuming any knowledge of them. In 
(a) and(b) we have the entanglement detection for the first 
state, when the measurements are considered in a different 
order at the program. In (c) and (d) we have the entangle- 
ment detection for the second state. In (e) and (f) we have 
the same type of analysis for the third state. See the main 
text for details. 
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